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p ,i Abstract 

, Casimir energy is a nonlocal effect; its magnitude cannot be deduced from heat kernel expan- 

^T) i sions, even those including the integrated boundary terms. On the other hand, it is known that the 

divergent terms in the regularized (but not yet renormalized) total vacuum energy are associated 
^ ■ with the heat kernel coefficients. Here a recent study of the relations among the eigenvalue density, 

■ the heat kernel, and the integral kernel of the operator e~*^^ is exploited to characterize this asso- 
. elation completely. Various previously isolated observations about the structure of the regularized 

CN ! energy emerge naturally. For over 20 years controversies have persisted stemming from the fact 

I that certain (presumably physically meaningful) terms in the renormalized vacuum energy density 

■ in the interior of a cavity become singular at the boundary and correlate to certain divergent terms 
r~| ' in the regularized total energy. The point of view of the present paper promises to help resolve 

these issues. 
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I. INTRODUCTION 



The aim of this paper is to bring order and completeness into a welter of observations 
in the physics literature about various contributions to quantum vacuum energy (before 
renormalization) and their relations to the heat-kernel expansion. The primary tool is the 
systematic theory of Riesz means of spectral densities and their relation to the asymptotic 
expansions of various integral kernels associated with the central partial differential operator 
of the field theory (e.g., |20|); the Casimir energy, regularized by an exponential ultraviolet 
cutoff, is identifiable with certain terms in these expansions jl8|. Both total energy and 
local energy density are considered; for the latter, we concentrate on the singular asymptotic 
behavior near a boundary, which was the su bjec t of considerable calculational attention and 
physical controversy two decades ago (e.g., [l^, HH), when the Casimir effect was of great 
interest as a model of quantum effects in cosmology (and in the bag model of hadrons). 
Interest in Casimir energy has recently resurged because of greatly improved experiments 
[l^HEHEIi^ and hints of technological application With it, the controversies have also 
reemerged (e.g., [i^ 3^ 3^)- (More complete bibliographical citations appear later in their 
proper contexts.) 

Section |n] sets up notation and summarizes the needed background information about 
the asymptotic expansions of integral kernels and the Riesz means of eigenvalue densities. 
Section IIIII concerns the energy density of a scalar field, with emphasis on the instructive 
example of parallel plates with Robin boundary conditions; to avoid confusion in the quan- 
tum theory later, it is important to note that the scalar field's energy contains, already 
at the classical level, a term concentrated on the boundary, unless the conf or mal- coupling 
parameter is given the nonstandard value |. The heart of the paper is Section HV] where the 
structures of the asymptotic expansions of the heat kernel and the total regularized energy 
(as functions of an auxiliary "time") and of the renormalized energy density (as function 
of the distance from the boundary) are presented and then correlated in dimensions 1, 2, 
and 3. In Section |V] we review the controversies over infinite boundary terms in the energy 
and draw some conclusions on the basis of the observations in the previous section. 



II. NOTATION AND FRAMEWORK 

Let H he a. self-adjoint, elliptic, second-order differential operator, independent of t. (To 
avoid some extraneous complications later, let us also assume that the spectrum of H is 
nonnegative.) Then the (hyperbolic) wave equation 

defines a field theory with time-translation invariance. Here we consider only scalar fields, 
but the extension to, for instance, the electromagnetic field would be straightforward. Also, 
we are primarily concerned here with cavities in fiat space, without potentials, so that H 
is just the negative of the Laplacian equipped with certain boundary conditions; however, 
most of what follows applies with only minor complications to more general operators 

H = -g(xrw^W^ + V{x) 

with gravitational, gauge, and scalar potentials. 
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For present purposes, we take the regularized vacuum energy of the quantum field to be 
defined as 

^W^^E^ne-"*, (2.2) 



n 



where Un are the eigenvalues of H, and t is a parameter with dimensions of time — not 
the physical time in ()2.H) . (Of course, at this point one is assuming that H has only 
discrete spectrum. Continuous spectrum is normally associated with an unbounded region 
in space, for which one would not expect a finite total energy to exist.) The physical 
(renormalized) energy is identified with the constant term in the asymptotic expansion of 
E{t) as t — ^ 0, in keeping with the expectation that each mode of the field with frequency 
u should contribute energy ^u. Formula ()2.2|) is a standard, but admittedly arbitrary, and 
recently somewhat unfashionable, way of regularizing the energy. The relation of this and 
other "ultraviolet" or "cutoff" regularizations to the now more popular "analytic" methods 
(notably dimensional reg ularization and zeta functions) has been clarified by several authors 



In quantum field theory one has not only a total energy but also a local energy density, 
Too(x) (and, indeed, an energy- momentum (stress) tensor, T^^{x)). The energy density 
remains meaningful in unbounded regions, where the total energy is not defined (and H may 
have continuous spectrum). There is a formal expression for Too (a;) as a sum or integral over 
the spectrum, the summand being built out of the eigenf unctions, or the spectral projection 
kernel, and their derivatives at x; for the details we refer to the literature (e.g., [3|) and 
Sec. mil From that expression a regularized quantity Too(a;;i) is defined by inserting an 
exponential cutoff in precise analogy with (j2.2p . In fact, it is known that this approach 
is not quite adequate when spatial curvature, or indeed any potential, is present, because 
trace anomalies arise; the conceptually closest adequate approach is then covariant point 
splitting (e.g., [6]). That complication is irrelevant, however, to cavities in fiat space without 
potentials, which are our concern at present. 

Let d be the spatial dimension, let M denote the d-dimensional region in question (where 
H acts), and let t be an auxiliary variable ranging from to +oo. The well known heat 
kernel, K{t, x, y) = (x|e~*''^|?/), solves the heat equation associated with H, in the sense that 

u{t,x) = / K{t,x,y)f{y)dy 

J M 

is the unique solution of the (parabolic) initial-value problem 

— = -Hu, u{0,x) = f{x). 
It is well known that K possesses an asymptotic expansion of the form 

oo 

K(t)~^6,r^+i (2.3) 

More precisely: (a) If a; is a point in the interior of M, then the diagonal value K{t, x, x) 
has an expansion of form (j2.3|) . where bs depends on x. (b) When M is bounded, the 
"trace" J^j K{t,x,x) dx has an expansion of form ()2.3|) . (The second statement is not a 
trivial consequence of the first, because expansion (a) is nonuniform in distance from the 
boundary. In fact, for a fiat-space cavity bs{x) in (a) is zero for all s > 0, but the same 
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is not true for the number bs in (b).) For a cavity, 60 in (b) is proportional to the volume 
of M, and the higher coefficients are integrals over the boundary of M whose details depend 
on the boundary conditions and on the curvature of the boundary. For a more general H, 
similar local integral formulas for bg in terms of the potentials, curvature, etc. inside M 
and on its boundary are known in great detail for small values of s, their calculation for 
ever-increasing s being a major mathematical industry (e.g., j32i] ) . 

The less well known cylinder kernel, T(t,x,y) = {x\e~^^\y), can be defined similarly: 



u{t,x) = / T{t,x,y)f{y)dy 
is the unique bounded solution of the (elliptic) boundary-value problem 

— = +Hu, u{0, x) = f{x). 

(The region in {t, a;)-space is a semiinfinite cylinder with base M.) The counterpart of ()2.H|1 
is 

00 00 
Tit)r^Y.^st-''^'+ E (2.4) 

s=0 s=d+l 
s—d odd 

(Again, ()2.4p has both (a)- and (b)-type interpretations.) The central fact relating these 
expansions (which in essence appears in [ll|. Sec. Ill]) is: 

Theorem. 

1. If s — is even or negative, 

e.--"=2-rg-| + i)., (2.5) 

(and = 0). 

2. If s — (i is odd and positive, 

= (_i)(-rf+i)/2^^i/22^-.+i^_l___ ^ (2.6) 

■'^ V2 ~ 2 2 

whereas Cg is undetermined by the heat kernel expansion. 



Before explaining where this theorem comes from, we make some remarks upon it. First, 
in ()2.5|) and ()2.6|) the gamma functions have no zeros or poles for the relevant values of 
the parameters; therefore, a zero or nonzero bs is always associated, respectively, with a 
zero or nonzero value for the coefficient on the left-hand side of the equation. Second, the 
"new" coefficients {s — d odd and positive) are inherently global in their dependence on 
the geometry of M: There are no local integral formulas like those for bg (and hence for 
the left-hand sides of ()2.5j) and ()2.(j|l ) fj^]. Third (see Sec. llVjl . the renormalized vacuum 
energy is equal to — |ei+rf (the first of the new, nonlocal coefficients) from expansion (b), plus 
possible local terms. (Note that a purely local formula for vacuum energy would imply that 
the energy associated with parallel plates is independent of their separation, hence that there 
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is no Casimir force!) Fourth, for a particular H any or all of the logarithmic coefficients fs 
may happen to vanish, but the nonlocal, independent nature of the corresponding Cg persists 
in such a case. 

The theorem can be proved by relating the two expansions, (j2.3j) and (j2.4|) . to certain 
"moments" (in a generalized sense of the word) characterizing the asymptotic behavior of the 
eigenvalue density (or of a local spectral density, the inverse Laplace transform of K{t, x, x), 
in the "untraced" (a) situation). One suitable set of such quantities consists of residues and 
values, at certain points, of the zeta functions associated with the operators H and ^/H. In 
that case the proof is an adaptation and generalization of the proof of Lemma 2.2 of |.2L] . 
(See also (HI.) We omit the details, citing only the crucial equation 

C^i2s) = Ch{s). (2.7) 

An alternative set, with more direct spectral and physical significance, consists of the 
Riesz means of the eigenvalue density. These are the rigorous counterpart of the formal 
high-frequency asymptotic expansions of the eigenvalue density whose term-by-term Laplace 
transforms yield the expansions ()2.H|1 and ()2.4|1 . The importance of Riesz means in under- 
standing ()2.4|) ■ and vice versa, has been developed in ISIl^l- (Papers [3, [3 113 
related.) The present paper demonstrates the usefulness of that work in physics. We shall 
review it here, u sing the language appropriate to the "traced" (b) situation. 

Riesz means |27l lisl 2^ generalize the Cesaro means used in place of partial sums to 



accelerate the convergence of Fourier series. Let /u(A) be the number of eigenvalues of H 
less than A (the counting function, or "spectral staircase"). Write 

X = uj^ (cu>0). (2.8) 

(Thus, if A is an eigenvalue, u is the frequency of the normal mode(s) of the wave equation 
associated with A; also, u is the corresponding eigenvalue of \/H.) For a positive integer a 
the Riesz mean i?^/i(A) is defined by averaging fi over an a-dimensional simplex: 

i?XA) = ^A-"|'---|MA)A (2.9) 

where a successive indefinite integrations are implied. (This operation can be reexpressed 
as a single integral with a nontrivial kernel function.) A different sequence of Riesz means 
is defined by integration over co: 

Ry^u;) ^ -u;-" j ... J ^{Co')du, (2.10) 

Starting from the known properties of the heat-kernel expansion (|2.3p . one can prove the 
finite asymptotic approximations 



a 

w — ^ ^ d—s „ / d — a — 1 \ , . 

i?> = E «a.A— + O (A^) , (2.11) 

s=0 

a 

Ry = y^Co,sOj'^-'+ V d^sUj'^-'\nuj + 0(uj'^-''-^\Yiuj) . (2.12) 



=0 s=d+l 
s—d odd 
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Here the a coefficients are related quite directly to the heat-kernel coefficients: 



d+s _|_ 1 
2 



Tis + 1) 



ass. (2.13) 



(The a^s with a ^ s have similar formulas, which we ignore here because those coefficients 
contain no additional information.) Similarly, the c and d coefficients are related to the 
cylinder expansion by [ip = (InF)'] 

Cg = — Css if (i — s is even or positive, (2-14) 

fs = - ^ ' = ^ [css + V^(c? + l)4s] if - s is odd and negative. 

The final link is the connection between (c^s, dgs) and a^s . In j20| the ratio of Css to 
was given as a complicated finite sum of gamma functions. However, hidden in the appendix 
of |2i| (because the authors did not appreciate their significance at the time) are formulas 
that, after some routine manipulation, yield the relatively simple result 



r(f-| + i)r(f + | + i 
2^r(f + |)r(f + i) 



ass if c? — s is even or positive. (2-16) 



When d — s is odd and negative (when the ffist factor in the numerator of (j2.16p has a pole), 
[ioj l shows that Css is undetermined by ass while dss comes into existence: 

(-i)d+i r(|-f) 

dss = ; rr-rr — ^ if (i — s is odd and negative. (2.17) 

(s-o?-l)!(i! r(-f-f) 

Equations (jTHl) and now follow from ()2.13|) - ()2.17|) with aid of 

r(^)r + = v^2i-2^ r(2z), r(-^)r(z + 1) = . (2.18) 



Note that the Riesz means ()2.9p and ()2.10j) provide two families of definitions of an 
"averaged density of states" with certain weightings. More precisely, or somehow 
describes the averaged asymptotic behavior of the eigenvalue density /i'(A) at high frequency, 
but it also contains constants of integration summarizing what happens to /x'(A) at low 
frequency. (As a increases, the smoothing out of the eigenvalue distribution becomes more 
drastic, so that the expansion p.llj) or ()2.12|1 can be meaningfully extended to a higher order; 
but at the same time, the contribution from the lower frequencies becomes more dominant.) 
The crucial fact is that the coefficients in the omega-mean expansion (equivalently, those in 
the cylinder-kernel expansion) contain more information than the coefficients in the lambda- 
mean (or the heat-kernel) expansion. Specifically, Css for s — d odd and positive is an 
independent, nonlocal spectral invariant. In particular, s~d = 1 corresponds (modulo local 
terms) to the Casimir vacuum energy. 
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III. ENERGY DENSITY, ROBIN BOUNDARY CONDITIONS, AND PARALLEL 
PLATES 



The energy density of a scalar field in fiat space-time is (e.g., [^) 



Too-- 



2 



f j +(l-4O(V0)^-4e0V^ 



(3.1) 



where ^ is an arbitrary number called the conformal coupling constant. In curved space- 
time ^ derives from a term (in the Lagrangian or the equation of motion) coupling (p to the 
curvature; in a gravitational context, therefore, different values of ^ correspond to different 
physical theories. In fiat space, however, ^ does not appear in the equation of motion, nor 
in the spectral decomposition of H, and the terms proportional to ^ in ()3.H1 form a total 
divergence, — 2^V ■ (0V0). The facile conclusion is that all values of C, yield the same total 
energy and hence are physically equivalent. We shall see, however, that subtleties arise at 
both classical and quantum levels. 

Because it leads to nontrivial phenomena even in one dimension, it is instructive to study 
the Robin boundary condition, 

l-.^. (3.2, 
Here 7(x) is, in general, a function defined on the boundary of the region M, and -k- is the 



outward normal derivative. Dimensionally, 7 is an inverse length. Romeo and Saharian |41 
have recently made an exhaustive study of the Casimir-like situation of two parallel plates 
with a constant value of 7 on each plate. (A continuing series of papers by Saharian and 
coworkers address similar issues in more complicated geometries). Of course, the special cases 
of the Neumann (7 = 0) and Dirichlet (formally, 7 = 00) boundary conditions have been 
understood for a long time. With the sign convention in ()3.2p . our condition of nonnegative 
spectrum for H is guaranteed by taking 7 negative; Romeo and Saharian consider both 
signs of 7 (which they call /3~^), ignoring the contributions to the energy from modes with 
negative eigenvalues, which do not affect the renormalization theory. 

If M is bounded, one should be able to integrate Too over M to obtain a conserved energy. 
(For parallel plates, integration in the perpendicular dimension yields an energy per unit 
plate area.) In fact, in the Robin theory it is necessary to include in the energy a surface 
term concentrated on the boundary: 



E = 



f Too(x) d^x + ^(4^ - 1) / 70' dS (3.3) 

JR4 2 JdM 



(where dS is the element of {d — l)-dimensional surface area on the boundary, dM). A 
calculation, using the equation of motion ()2.1|1 for cj) and an integration by parts, shows that 
^ = 0, the boundary term being essential for this result if 7 7^ and ^ 7^ |. Such a term 
was apparently introduced into quantum field theory — through a different rationale — in 
[3]| . Classically, the term is presumably related to the fact that a one-dimensional Robin 
problem models a vibrating string attached to a point mass, which can exchange energy 
with the string (e.g., [l3, p. 15]). 

Three things should be noted about the boundary term in ()3.3p . First, it is a completely 
classical matter; in particular, it is not a renormalization counterterm (at least, not yet). 
Second, it arises only in the full Robin theory: it vanishes in the pure Neumann case (7 = 0) 
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and the Dirichlet case (0 = 0, where the conservation of E without that term can be easily 
verified). Third, even in the Robin case it vanishes if ^ = |. Changing ^ (in fiat space) 
amounts to redefining some of the boundary energy as spread throughout the interior of M 
by integration by parts, and | is the choice that moves the energy entirely into the interior. 
The striking thing about this value of ^ is that it is neither the "minimal" choice [C, = 0) 
nor the choice that makes the full theory, including gravitation, conformally invariant. The 
"conformal" value is well known to be 

which equals when d = 1, equals | when d = 3, and approaches | as d — >■ oo. 

In the quantum field theory one evaluates the expectation values of observables such as 
Too{x) and E in the vacuum, or ground state, of the field. The general rule of thumb 
for doing so, when the observable is classically defined by a quadratic form T(0, 0), is to 
evaluate the associated bilinear form on the normal modes of the field and sum: 

{0\T\0)=J2Ti4>nA*n)- (3.5) 

n 

Here (f>{t,x) = (2ci;„)^^/^0„(a;)e~*'^"*, where 0„ is an eigenvector of H with eigenvalue Un] 
generalization to continuous spectrum is possible. (Henceforth the vacuum bra-ket symbols 
will be omitted when there is no danger of confusion.) The result of this formal calculation 
is usually a divergent sum or integral, which must be regularized — perhaps by inserting 
an exponential ultraviolet cutoff, exactly as in ()2.2j) — and ultimately renormalized. (How- 
ever, the calculations of energy density in |^ do not require an explicit regularization step; 
instead, the renormalization is accomplished by isolating and subtracting at the integrand 
stage the terms corresponding to the vacuum state in infinite fiat space.) When such cal- 
culations were done some 20-30 years ago for the Dirichlet and Neumann problems with 
parallel plates, it was found that when ^ = the (renormalized) energy density Tqq{x) is 
constant (and nonzero) everywhere between the plates, thus yielding a finite Casimir energy 
(per unit area) that must be attributed to the global configuration of the plates, not to an 
interaction between the field and the boundary per se. (In the Neumann case the contribu- 
tion of the mode with a; = is being ignored here. When a zero-frequency mode exists, there 
is no true vacuum state, but vacuum conditions can be approached arbitrarily closely.) On 
the other hand, for ^ 7^ the renormalized energy density includes a contribution strongly 
concentrated near the plates themselves; that part is, in fact, nonintegrable, hence at least 
superficially inconsistent with a finite value for the renormalized total energy. When the 
boundary is curved, there are such near-boundary terms even when the coupling constant 
takes the conformal value, although the singularity at the boundary is weaker in that case. 
The recent Robin results j4l| are a consistent extension of that picture: 

• For a single (fiat) plate, the renormalized Tqq vanishes if ^ = ^c, regardless of 7. 

• In general, the energy density around a single plate is proportional to ^ — ^c- If a; is 
distance from the plate, Tqq(x) is proportional to x"^'^"''^^ in the Dirichlet and Neumann 
cases, and in the general Robin case it has additional terms proportional to ^^^^'^"'^'^"'"^^ 
{k = 1, . . . ,d) and In |7x|. (Romeo and Saharian actually find an exact expression for 
Too, of which these terms are merely the divergent part of the asymptotic expansion 
for small x.) 
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• For two parallel plates, Too is constant in the space between the plates if ^ = ^c, with 
a rather complicated dependence on the parameters 071 and 072 (a being the plate 
separation and 7^ the Robin coefficients on the two plates), which reduces for the pure 
Dirichlet or Neumann problem to the known a~^^~^^'^ dependence. 

• For parallel plates with ^ 7^ ^c, the result is even more complicated, but it can be 
written as the sum of the effects of the individual plates (described above) plus an 
interaction term that is everywhere finite. 

• Romeo and Saharian perform an independent calculation of the renormalized total 
energy by the zeta-function method. This calculation is independent of ^, as the 
eigenvalues do not depend on ^. The result includes a term that depends logarith- 
mically on an arbitrary scale parameter — indicating that a logarithmic divergence 
would appear in some other regularization methods. This term is independent of a, 
hence presumably unobservable. 

• The authors demonstrate a certain formal consistency between their energy-density 
and total-energy calculations. In particular, for a single plate they find 

^vol = / Too(x) dx = 4t/(e - ^c)^, ^surf = "(4^ ' 1)7 (0 1 0(0)^ 1 0) = - 40^, 
Jo 2 

(3.6) 

where E = E^oi + -E'surf is the putative total energy per unit area (on one side) of the 
plate. However, they point out that the integrals defining the three quantities in ()3.6|) 
are actually divergent. One can eliminate the volume term by taking ^ = , or one 
can eliminate the surface term by taking ^ = | , but not both at once. 

In view of the last point, a strong case can be made that the most convenient choice of 
^ is ^ , rather than the traditional values, or ■ One can then avoid accounting for and 
renormalizing the surface term, and one can be confident that any volume contribution must 
be taken seriously, not dismissed as an artifact. 

We have summarized the paper jiJ] at some length because the parallel-plate Robin 
system is such an excellent, explicitly solvable model of the phenomena and issues that arise 
in more general systems, notably cavities with curved boundaries. The principal extrinsic 
curvatures of a boundary surface have the dimension of inverse length and enter calculations 
or expansions of energy density and total energy in very much the same way as does the Robin 
coefficient 7. The leading divergence in Too as the boundary is approached is ubiquitous when 
^ does not have its conformal value; unlike the Dirichlet case, the classical surface term in 
the total energy of the Robin system interferes with the facile dismissal of this problem 
as some kind of artifact of a poor definition of energy density (whereas the nonexistence 
of the surface term in the Dirichlet case shows that the problem can't be resolved by a 
simple cancellation of volume and surface contributions, either). The higher-order terms in 
the expansion of Too near the boundary, which appear for any ^ when 7 7^ 0, are closely 
analogous to terms that appear near curved boundaries for the electromagnetic field, as well 
as the conformally coupled scalar one. 

IV. SCHEMATICS OF THE EXPANSIONS 

Consider now a general cavity whose boundary is characterized by an extrinsic curvature 
tensor, k{x), and also equipped with a Robin function, 7(0;). The structure of the heat 
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expansion ()2.3p . for the trace, is 

K{t) ~ r^'^iV + St^l^ + (k + 7)t + (/t' + 7' + '^7)^'/' + + O^' + ■■■]■ (4.1) 

Here, and in all similar equations later, all the terms are indicated very schematically: all 
numerical coefficients, the integrations over dM, and the tensorial nature of k when d > 3 
are suppressed, as are terms involving derivatives of 7 or k (which are not important in 
low orders and dimensions). V is the volume of M, and S is the surface area of M (more 
precisely, the {d — l)-dimensional volume of DM). 

The schematic form of the (traced) cylinder expansion ()2.4|) now follows by the theorem 
in Sec. IHI The significance of the structure depends on the dimension. In the following 
formulas, the coefficients in boldface are the new, nonlocal terms, , and those in square 
brackets are constant terms, which do not contribute to the regularized energy, 

, 1 dT 

m = -^^. (4.2) 

(The trace of the operator e~*^ is Y.n e"'^"*, so ()4.2p is equivalent to ()2.2j).) 

d=l: T{t)r^Vt-'^ + [S]+-ft\iit + Et + -fH'^ + -fhHnt + Ft^ + --- . (4.3) 

(In dimension 1 there is no curvature, and the "surface area" S degenerates to the number 
of endpoints of the region — usually 2. The 7 terms are also sums over the endpoints.) 

d = 2: T{t) ~ Vt-'^ + St-^ + [K + -f] + {K'^+-f'^ + K-f)t\nt + 'Et + {K^ + ---)t'^ + --- . (4.4) 

d = 3: T{t) ^Vt-^ + St-'^ + {K + -i)t-^ + [K^ + -i'^ + K-i] + {K^ + ---)t\nt + 'Et + --- . (4.5) 

The untraced kernel expansions can be discussed similarly, but there is less to say. In the 
cavity setting, where there are no potentials or spatial curvature, the local heat expansion 
is trivial: 

i^(t,x,x) ~r'^/2^0(t°°). (4.6) 
The local cylinder expansion, therefore, has the structure 

T(t,x,x) ~t""' + E(x)t + --- . (4.7) 

Let Too(t, x) be the energy density, regularized by an exponential cutoff as discussed below 
()3.5|) . The relation between Too(t, x) and T(t, x, x) is less straightforward than that between 
their global counterparts, E{t) and T[t). The analogue of ()4.2|) yields only the terms in ()3.ip 
involving time derivatives. The other terms can be obtained [l8| by applying appropriate 
spatial partial derivatives to the integral kernel of the operator 

^ e-*^, (4.8) 



H 

which has not been as extensively studied as the cylinder operator but has a very similar 
nature. (Note that the original cylinder kernel is the time derivative of this new one.) For 
the qualitative considerations of the present paper it is not necessary to delve further into the 
details; it is not misleading to think of Too(t, x) as "essentially" the t-derivative of T(t, x, x). 

This is as far as the heat kernel can take us. To find E(x) one must do independent 
calculations. But the literature, notably [13, i^l, provides considerable information 
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about the behavior of the renormahzed Too (a;) as x approaches a boundary. Our task now 
is to compare these expansions in distance from the boundary (which we once again denote 
by x) to the expansions (in t) of the regularized total energy, the derivatives of ()4.3|) - fj4.5p . 
The point is that the finite, renormalized energy density generally is not integrable up to 
the boundary, and if such divergent terms are taken seriously (and not compensated by 
postulated terms concentrated exactly on the boundary j3l|, whose discussion we postpone 
to Sec. 13); then for consistency they must somehow be reflected in terms of E{t) that become 
infinite in the limit t — *• and hence would be removed in a conventional renormalization of 
the energy. 

A. Dimension 1 

The energy density found by Romeo and Saharian jij] near a Robin endpoint (with 
^ = i) can be correlated with the derivative of ()4.3|) as follows: 

S 7 

Too ~oo + — + — +7^1nx+ finite 

i i i \ i X (4.9) 

The leading oo in Too is the universal, x-independent formal vacuum energy of infinite empty 
flat space, which everyone agrees should be discarded in renormalization. It is included 
here because it clearly corresponds directly to the volume divergence in E. The arrows 
indicate how other terms in the two series might be related, on the basis of their dependence 
on 7. (Whether all these arrows represent genuinely existing relationships remains to be 
investigated by a thorough and rigorous calculation of Too(t, x, x), followed by an integration 
over X with t finite but approaching 0. That project is beyond our present scope.) The most 
striking observation is that the "surface" term in the density has completely disappeared 
from the total energy. (Recall that the corresponding term in ()4.3|) was destroyed by the 
differentiation.) The 1/x term in the density might be held responsible for the \nt term in 
the energy, as well as part of the finite term. The Inx term is integrable, and its integral 
is presumably included in E. Finally, we note that all the explicitly 7-dependent terms 
in E will be independent of the length of the cavity, so they will not contribute to the 
Casimir force between the endpoints (usually considered the only quantity of experimental 
relevance). The force must come entirely from E. 

B. Dimension 2 

From now on, for brevity, the Robin function, 7, will not be explicitly indicated in the 
displays (j4.1Up and (j4.1ip . because it enters in the same way as the curvature, k. I do 
not know a reference for the dependence of Too(a;) upon k in dimension 2, but it appears 
to be dictated by dimensional analysis and analogies with the Robin formula j4l| and the 
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three-dimensional results 



Too ~ oo + — + — + — + K^lnx + finite 

i i i i \ i X (4.10) 

V s 

E T + ^ + K^hit + + k'^) + kH + ■■■ 

Most noteworthy here are the term k'^/x in the density and the corresponding term K^lnt 
in the total energy. The analogous terms in dimension 3 are proportional to k^. When 
total energy calculations are done for both the inside and the outside of a boundary surface 
(a "shell" modeling a thin electrical conductor), it has often been noted that divergent 
contributions proportional to odd powers of k, will cancel, because k has opposite signs 
on the two sides. This fact gives rise to an important difference between even and odd 
dimensions. In even dimensions, the Int term in E is proportional to an even power of 
K, SO it survives even for thin shells. Logarithmic terms in E{t) regularized by ultraviolet 
cutoff are precisely the divergences that persist in analytic regularization schemes, giving 
rise, under even the most liberal interpretation 0|, at least to ambiguous finite terms in the 
renormalized energy that depend logarithmically on an arbitrary mass scale. In dimensional 
regularization, the k'^ divergence shows up as a pole at each even value of d in the total 
Casimir energy associated with a spherical shell [35|. Here we have seen in a more direct, 
elementary way that such a term is dictated by the known form of the heat kernel. 



C. Dimension 3 

We are now in the domain of the classic results of Deutsch and Candelas and Kennedy, 
Critchley, and Dowker [31], which also apply to the electromagnetic field. 

Too~oo + — + — + — + — + O(lnx) 

OC 3j 3j 

i i i I i \ i (4.11) 

E - J, + J, + J2 + ^1^^ + + + ■■• 

Let us review the litany of special features that make the Casimir force on a spherical 

conducting shell finite, even when calculated by integrating the energy density right up 
to the boundary: The volume term is removed by the local renormalization. The surface 
term vanishes for electromagnetism because of a cancellation between electric and magnetic 
contributions. The k and tc" terms cancel (between inside and outside) for any thin shell. 
That leaves . In dimension 2, k is a symmetric matrix with two eigenvalues, ki and K2 (the 
principal curvatures) . There are thus two independent quadratic invariants to be integrated 
over (9M, which may be chosen as 

D = KiK2, /\={Ki-K2f. (4.12) 

For a sphere, A is identically zero. For any closed surface, Jq^j D{x) dS is a topological 
invariant by the Gauss-Bonnet theorem, so its derivative with respect to the sphere's radius 
vanishes. (More generally, for any bounded region M, any quadratic integral of the type 
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IdAt is constant under a shape-preserving magnification of M, even if A is involved. 

Note that this observation does not settle how the energy per unit length of an infinitely long 
cylinder depends on the radius, since the scaling of the area is different in that case.) As 
emphasized by Candelas 0, H, [13 > one would nevertheless expect infinite energy differences 
to result from shape-changing deformations of a closed surface, because of the A term. It 
is therefore significant that no term of the type n'^/t can appear in the regularized total 
energy, E{t), because of the structure of the traced cylinder expansion, ()2.4|1 or ()4.5|1 . We 
return to this issue in the next section. 



V. IMPLICATIONS 

The theory of vacuum energy in electromagnetism has increasing experimental verification 
and many indirect implications for observable phenomena. Nevertheless, the simplest and 
crispest theoretical models tend to be idealizations that are not directly testable. Perhaps 
for this reason, controversies persist about renormalization procedures and the physical 
significance of various divergences. These disputes fall into two main categories. 

First, there is the broad question of whether "logarithmic" divergences are a sign of a bad 
model, or just an inherent renormalization ambiguity. Among the divergences that arise in 
energies regularized by point-splitting or ultraviolet cutoffs, such as the second line of (j4.11|) . 
only those proportional to In t survive 0, El in analytic regularization schemes, where 
they appear as poles at the physical values of the regularization parameters. The analytic 
methods therefore give finite energies automatically in many problems, without a need for 
a renormalization subtraction, but in some other problems are embarrassed by a surviving 
divergence whose removal introduces an unavoidable ambiguity parametrized by an arbitrary 
length or mass scale. Some authors 13. l3]| regard this as not at all a problem, while others 
regard the poles as "truly disturbing" 'Sa', p. 199] or at least puzzling. 

Second, and specific to Casimir calculations, is the question of whether the nonintegrable 
divergences in the renormalized energy density near boundaries must be taken seriously, even 
in cases when the renormalized total energy (or the observable force) is finite. In fact, how 
can such a divergence even be regarded as logically consistent with a finite total energy? 
Both Kennedy et al. [sH and Milton and coworkers [s^, 113] trust the calculations of 
finite total energies and forces more than the calculations of local energy densities. Kennedy 
et al. argued that the infinities in Tqo{x) must be compensated by delta-function terms 
concentrated exactly on the boundary itself, so that the finite total energy predicted by zeta- 
function calculations will result. (The boundary term in ()3.3|) originated in this discussion 
in [31], without a clear explanation of how it could be relevant when ^ = 1 or 7 = or 00.) 
Milton writes j3^, p. 237], ". . . [T]he local vacuum-fluctuation energy density is to a large 
extent meaningless.. . . [D]ivergences in the energy which go like a power of the cutoff are 
probably unobservable, being subsumed in theproperties of matter." 

On the other side, Candelas and Deutsch p, H, [13 argued that the boundary infinities 
in the energy density are real (within the theory) and must be regarded as a defect of 
an overly idealized model. T his p osition has recently been reexpressed by Graham, Jaffe, 
and coworkers ji^ [2^ [2^, [2^ [3^. The argument, stated in electromagnetic terms, is that 
no real material is perfectly conducting at arbitrarily high frequencies. When an idealized 
boundary is replaced by a more realistic model of the interaction of the quantized field 
with other matter, the infinite pileup of energy near the boundary will become finite, but 
possibly very large and definitely physically real. A strong argument for this position is that 
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the energy-momentum tensor acts as the source of the gravitational field, so a nonintegrable 
singularity in it is physically unacceptable, even if the total energy comes out finite by virtue 
of a cancellation (as happens in some scenarios involving thin shells). 

The present paper does not resolve all these issues, but it casts some light on them and 
suggests some directions for future research that will clarify them further, as indicated in 
the following remarks. 

A. Thin shells; Robin divergences 

Much has been made of the disappearance of some divergences by cancellation in the 
case of thin shell boundaries. The dimensional dependence (even vs. odd) of this effect is 
quite clear from (j4.9|) - (|4.11|) and their antecedents. Those formulas also indicate analogous 
divergences specific to boundaries with nontrivial Robin coefficients, which were discovered 
in j4l| but could have been anticipated from an understanding of the relation between 
vacuum energy and the heat kernel. Note that there is no reason to expect any cancellation 
of exterior and interior Robin effects for a thin boundary. 

B. Disappearing terms; total derivatives; conductors with anisotropic curva- 
ture 

A pervasive feature of ()4.9p - ()4.1H) is that terms in Too proportional to {x being 
distance to the boundary) never have counterparts in the regularized total energy. On 
grounds of dimensional analysis such terms would be expected to appear in E{t) and be 
proportional to t~^, but the derivation of E{t) from the heat and cylinder kernels shows 
that they are not there! (The crucial fact is that terms proportional to t^lnt occur in the 
cylinder kernel only if n is odd and positive — in particular, not for n = 0.) This absence of 
terms was noted by Cognola et al. and very recently by Bernasconi et al. 0, who 
found that it persists when the exponential cutoff is replaced by a Gaussian one. 

In dimension 1 the term is precisely the one associated with nonconformal coupling 
— that is, the leading boundary divergence, which (in any dimension, but for a scalar field) 
does not occur for ^ = and arises for other values of ^ from a term proportional to 
(^ — ^c)V ■ (0V0) in the formal energy density (cf. (j3.1|) and ensuing discussion). Since this 
term is a total derivative, it ought to give when integrated. Effectively this seems to be true, 
since the renormalized total energy is finite and independent of ^. As has been recognized (if 
not "understood") since jl^ and [s^], integration over x somehow does not commute with 
renormalization. There is a simpler and completely understood analogous situation with 
the heat kernel: The local heat expansion ()4.6|) is nonuniform in x, and hence integrating 
it over M does not yield the correct trace expansion, ()4.H) . The interesting terms in the 
latter come from a part of K{t, x, x) that concentrates entirely on the boundary in the limit 
t ^ 0. It is tempting to think of the boundary terms in Too (see (l4.9j) - (j4.1H) ) as the cylinder- 
kernel analogues of these boundary terms in K, but they are not: they are nonvanishing 
even when t = 0. More likely, there is some other part of T{t,x,x) that concentrates right 
on the boundary in the limit and partly or wholly compensates the integrals of the energy 
divergences near the boundary. This would corroborate, within the framework of ultraviolet 
regularization, the proposal of Kennedy et al. |3l|, made in the framework of zeta-function 
regularization. (It would not, however, solve the problem of a singular Einstein equation.) 
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We propose, therefore, a program to calculate T{t, x, x) in various models for finite t 
and to integrate it over x, and then to examine the limit t ^ 0. The expectation is 
that the results (including contributions from the kernel of fj4.8|) if necessary) will be fully 
consistent with the direct regularization and renormalization of the total energy, and that 
the details of the process will demonstrate precisely how the terms in the expansion of the 
renormalized Toq{x) near the boundary are related, consistently, with the divergent terms in 
E(t). Some steps in this program can already be extracted from the literature: Parentani 
j40| has made a similar observation in the context of a uniformly accelerating boundary. 
Ford and Svaiter also with an eye to resolving the paradox of the total energy of 
the nonconformally coupled scalar field, studied a statistically fluctuating boundary. Their 
formulas seem mathematically identical to the results of ultraviolet regularization with an 
arbitrary (not necessarily exponential) cutoff function, with a fluctuation parameter in the 
place of t. They find (in dimension 3) that the nonconformal term of Tqq oscillates from 
negative to positive and back to negative, in such a way its integral is indeed zero. Finally, 



Graham et al. [23|, |2J, |26|, |39[ replace the ideal boundary with an extended object modeled 
by a second quantized field. Their results are very similar to those of |16| in dimension 3 and 
have one fewer oscillation in dimension 2. (Olum and Graham j3^ also point out that the 
boundary itself, as a physical object, has an energy density of its own that can be expected 
to swamp the energy induced in the scalar field; nevertheless, they conclude that the total 
energy density remains negative in a small region displaced slightly from the boundary.) 
All the available results suggest that the cancellation is insensitive to the functional form 
of the ultraviolet cutoff, and also to its physical interpretation (fluctuating boundary, soft 
boundary, or mere mathematical device). 

In dimensions greater than 1 the leading boundary divergence associated with a noncon- 
formal stress tensor does not go as so the additional boundary terms that render the 
regularized energy ^-independent are not related (at least, not directly) to the x~^-versus-t~^ 
phenomenon. Instead, that phenomenon would appear to have some other implications. In 
particular, in dimension 3 it indicates that E contains no term proportional to A = (ki—k^)^. 
This conclusion seems to be robust against changes in the form of the ultraviolet cutoff j^l ■ 
But local density calculations indicated to Candelas and Deutsch 0, H, E^l that that term 
would be not only large but negative (persisting into the renormalized energy for a physical 
conductor), thereby rendering a conducting foil unstable against wrinkling. This may well 
be true for a physical, imperfect conductor. However, in the framework of formal ultraviolet 
regularization of the idealized theory, our conjecture is that there is a compensating posi- 
tive term that concentrates right on the boundary as the regularization parameter is taken 
to 0, thereby removing the apparent conflict between local and global calculations. Clearly, 
detailed calculations are needed to verify this conjecture and to establish its relevance to 
realistic conductors. (Incidentally, the sharp-eyed modern reader may notice that the orig- 
inal argument concerning the A term in Appendix A of 2| incorrectly assumed that it 
is correlated with a term of order uj~^ in a formal expansion of the counting function fi. 
Indeed, Candelas pointed out in a later paper 0, Sec. 7.3] that no such term exists in /x, but 
that its alleged effects can still arise in the integrated contribution of the oscillatory part 
of the eigenvalue distribution. In the terminology of our Sec. |nj the relevant terms in the 
cylinder kernel and in Too are associated with constants of integration in Riesz means rather 
than with the local behavior of /i'(A), or even /i(A), at large A. Therefore, this tangential 
technical issue in no way invalidates the physical conclusions of 
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C. Genuine surface effects 



On the other hand, almost certainly there are some real near-boundary energies. In ()4.9|) ~ 
()4.11|) one observes terms of order in the renormalized energy density, which appear to 
be related to the divergent terms of order \nt in the regularized total energy (although 
explicit calculations, of the sort described above, should be performed to confirm this). As 
previously mentioned, terms of the latter type are agreed to yield terms with inherently 
ambiguous coefficients in the renormalized energy. Also, terms of order In X in Too s-^s 
integrable and presumably contribute part of the finite term in E. 

D. Nonrelativistic Casimir energy? 

On p. 30 of [3^, Milton raises the question of why, mathematically, there is no Casimir 
energy in nonrelativistic quantum theory. (Physically, such energy can't exist because the 
relativistic effect vanishes exponentially with large field mass.) He mimics the relativistic 
zeta-function calculation — for parallel plates with Dirichlet conditions — and finds an 
expression containing a gamma function in the denominator, whose pole forces the quantity 
to vanish in the relevant limit. The same conclusion can be reached more transparently from 
our point of view. The formal energy is 



But the sum here is simply T{t), which has, as in (|4.9|) - (j4.1H) . the schematic form 

T{t) ~ divergent terms + Et + {k'^+^ H h 7^^+^)^^ + Ft^ H , 

the coefficient of the term being dictated by the heat-kernel expansion. Therefore, in the 
renormalized energy one gets 



because k and 7 are zero in the system studied. 

In summary, what we have shown is that various results for specific models that previously 
emerged from forests of gamma and Riemann zeta functions now are easy consequences 
of the known general forms of the heat and cylinder expansions. There are, to be sure, 
enough gamma functions in and Sec. [TTl to satisfy anyone, but they are handled once 
and for all in the general theory, rather than arising in each concrete calculation out of 
some physically opaque analytic continuation. With respect to the controversies over global 
surface-energy renormalization versus local energy- density calculations, we have not taken 
a position squarely on either side. Rather, we have presented evidence that each side is 
correct in certain respects. In particular, a mechanism of partial surface cancellation has 
been pointed out that shows promise of reconciling the two points of view. For the leading 




2m ' 



which can be regularized in analogy with ()4.2|) as 

— — Ve 

4m df^ V ^ 
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boundary divergence of the scalar field, this mechanism has been fairly well established by 
other researchers for the X ^ divergence, it is conjectured by the present author 

as the only visible way of restoring consistency between the mathematical structures of the 
local and global asymptotic expansions. 
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